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The fundamental question of how Bose-Einstein conden- 
sates tunnel into a barrier is addressed. The cubic nonlinear 
Schrodinger equation with a finite square well potential, which 
models a Bose-Einstein condensate in a quasi-one-dimensional 
trap of finite depth, is solved for the complete set of localized 
and partially localized stationary states, which the former 
evolve into when the nonlinearity is increased. An immediate 
application of these different solution types is tunable tunnel- 
ing. Magnetically tunable Feshbach resonances can change 
the scattering length of certain Bose-condensed atoms, such 
as 85 Rb, by several orders of magnitude, including the sign, 
and thereby also change the mean field nonlinearity term 
of the equation and the tunneling of the wavefunction. We 
find both linear-type localized solutions and uniquely nonlin- 
ear partially localized solutions where the tails of the wave- 
function become nonzero at infinity when the nonlinearity 
increases. The tunneling of the wavefunction into the non- 
classical regime and thus its localization therefore becomes 
an external experimentally controllable parameter. 



I. INTRODUCTION 

One of the most attractive aspects of trapped dilute 
gas Bose-Einstein condensates (BEC's) is the number of 
experimental parameters which can be varied over a wide 
range. For example, the effective dimensionality of a con- 
densate can be changed by varying the confining frequen- 
cies of a harmonic trapping potential [|l]j^]. The size of 
nonlinear collective excitations such as solitons ||,|| and 
vortices @|| , which is characterized by the healing length 
of the condensate , can be varied over many orders of 
magnitude by changing the density || or using a magnet- 
ically tunable Feshbach resonance to change the effective 
atomic interaction strength ||. The mean field of the 
BEC is characterized by a macroscopic wavefunction. ft 
is here demonstrated that, in a trap of finite depth, tun- 
neling and localization of this mean field wavefunction is 
an externally controllable experimental parameter. 

A highly successful mean field model for BEC's is 
the Gross-Pitaevskii equation ]Tc|Jll| |. In its quasi-one- 
dimensional regime, which holds when the transverse di- 
mensions of the condensate are on the order of the con- 
densate's healing length and its longitudinal dimension 
is much longer than its transverse ones Jl2|,[l3| and is 
the regime considered in this paper, the Gross-Pitaevskii 



equation is identical to the well-known cubic nonlinear 
Schrodinger equation (NLS). Conversely, and not consid- 
ered here, when the condensate's transverse dimensions 
are much less than the healing length then the Gross- 
Pitaevskii equation no longer applies and other physical 
models are required 

The NLS models many other natural phenomena as 
well, including spin waves in magnetic materials [ jl5| , 
Bose-condensed photons 16 1, disordered media jt7|, he- 
lical excitations of a vortex line |Q , and light pulses in 
optical fibers |jl9). However, as BEC's are trapped, an ex- 
tra potential term is required. This has opened up new 
theoretical and mathematical investigations |]l^ , po| -p3| . 
The stationary cubic NLS with a potential may be writ- 
ten n 

[-dl + rj | f(x) | 2 +V tra *> (x) } f(x) = n f(x) , (1) 

where f(x) is the longitudinal portion of the mean- field 
condensate wavefunction, fi is the eigenvalue and the 
chemical potential, and 77 is the strength of the nonlinear- 
ity, an extra parameter from the usual NLS which cannot 
be scaled out. All quantities in Eq. ([[]) are dimension- 
less. Physical units relevant to present experiments are 
r/ on aN, where a is the s-wave 



discussed in Sec. IV 



scattering length which characterizes binary atomic in- 
teractions in a dilute gas and N is the number of atoms 
in the trap. Note that the transverse degrees of freedom 
have been frozen out. 

In order to consider a BEC in a finite trap the square 
well potential 



ytrap^ 



\x\<l 

V \x\ > I 



(2) 



is studied. This potential gives analytic solutions, un- 
like harmonic traps, which in one dimension do not give 
rise to analytic solutions. It has the advantage of hav- 
ing a direct analog in the linear Schodinger equation, 
for which the stationary states have been worked out 
completely. In addition it models experimental config- 
urations which could be constructed with present BEC 
technology. 87 Rb BEC's trapped in hollow blue-detuned 
laser beams (2^] can be made quasi-one-dimensional [||. 
Laser light sheets already in use as endcaps for such cylin- 
drical traps have a variable intensity I which could make 
the trap finite. 85 Rb, for which the interaction length is 
particularly responsive to magnetic tuning via a Feshbach 
resonance p5[ , could also be trapped in this manner p6[ . 
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Solutions of Eq. (Q) with the finite square well potential 
Eq. (|^) are presented. When r\ is negative the solutions 
are localized. When r\ is positive and the chemical po- 
tential /i is less than or equal to the well depth V , the 
solutions are again localized. However, when rj is positive 
and /x > V , the solutions are partially localized, i.e. the 
tails of the wavefunction approach a constant non-zero 
value at infinity. These latter solutions grow analytically 
and continuously out of the localized, bound solutions, 
and have no analog in the linear Schrodinger equation. 

The transition between a localized and a partially lo- 
calized stationary state is an adiabatic one which does 
not change the number of nodes, unlike in the linear case, 
where the wavefunction aquires nodes when pushed out 
the top of the well. The nonlinear analog of scattering 
states, which continue to oscillate as x — > ±00 and do 
not grow analytically and continuously out of the bound 
state solutions, are not considered here. In the linear 
Schrodinger equation the square well parameters V a and 
I completely determine the stationary states, whereas in 
the NLS r\ is an extra external control parameter which 
determines the extent of the tunneling and therefore the 
localization of the wavefunction. Exploration of this ex- 
tra freedom is the key advance in this paper. 

The mathematical question of how solutions to the 
NLS decay under a potential step is, to the authors' 
knowledge, here answered for the first time. The clos- 
est previous treatments have considered a repulsive delta 
function, which models a small scatterer or impurity [ 27 , 
or approximate solutions to a harmonic potential [|28| 30 1 . 
Limits of the finite square well include the delta function 
and the infinite square well, the latter of which has been 
solved elsewhere [ |l2| , p0| ]. Both attractive and repulsive 
BEC's are treated, i.e. i] negative and positive. The 
problem of double well tunneling, which is well known, 
in the case of a high barrier, to lead to breakdown of the 
NLS |3l|]32|l , and therefore represents an entirely different 
physical situation, is not considered. 

The article is outlined as follows. In Sec. || the full set 
of localized and partially localized stationary solutions 
of the stationary NLS with a finite square well potential 
are presented. In Sec. Ill the delta function and infinite 
square well limits are explicated. In Sec. IV the appli- 



cation of these solutions to tunable tunneling in BEC's 
is discussed, with specific parameters relevant to current 
BEC experiments. Finally in Sec. [v| the results are briefly 
summarized. 



II. FINITE SQUARE WELL 

Just as for the linear Schrodinger equation with a fi- 
nite square well potential, solutions to Eq. (Q) with the 
potential Eq. (Q) are constructed by matching decaying 
solutions outside the well to periodic solutions inside the 
well, subject to continuity of f(x) and its first derivative 
at the edges of the well. Periodic solutions to Eq. (fil) 



with V(x) — are constructed from the special functions 
known as Jacobian elliptic functions |]33" J3~4| . Limits of 
these functions also constitute decaying solutions. Ellip- 
tic functions are more commonly known as solutions to 
the anharmonic classical oscillator, i.e. 9 + 9 — # 3 /3! = 0. 
The standard notation sn(a: | to) will be used for elliptic 
functions, where m is the elliptic parameter. This param- 
eter interpolates the Jacobian elliptic functions between 
trigonometric and hyperbolic functions. Their quarter 
period is given by the complete elliptic integral K(m). 
A short summary of the particular functions relevant to 
this work is given in Table |. 



) with V(x) = are 
12[. Of these, six are 



Trivial phase solutions of Eq. ( 
the set of twelve elliptic functions 
bounded and six are not. Of the bounded functions, three 
have different physical form, i.e. differ by more than a 
translational shift or a rescaling of the amplitude. These 
are cn(x | to), dn(a; | to), and sn(a; | to), and are all 
candidates for solutions inside the well. The to = 1 hy- 
perbolic limits of the remaining six unbounded functions 
constitute the tails outside the well. These are sech(x), 
csch(x), and coth(a;). to — 1 corresponds to the infinite 
wavelength limit of the elliptic functions. However, only 
certain combinations of these functions are possible. For 
i] > 0, sn(x | to) is the only periodic solution in the well 
and coth(a;) and csch(x) are both possibilites outside the 
well. For r\ < 0, cn(x | to) and dn(x | to) are both 
possible inside the well while sech(x) is the only possible 
decaying solution outside the well. Therefore there are 
four possibilities, two for each sign of rj. Of these four, 
two shall be shown to be localized states, one a partially 
localized state, and one not a solution. Thus there are 
three solution types. 

All solutions will be written in the form 



fi(x) x < — I 

fix) = { h{ X ) \x\ < 1 

f 3 (x) x > I 



(3) 



in correspondence with the form of the potential, Eq. (||) . 
These will be subject to continuity of f(x) and f'{x) at 
x = ±1 and, for the localized states, the normalization 
dx I f(x) 1 2 = 1. For the partially localized states, 
which are not normalizable, the normalization condition 
will be replaced with the use of the constant amplitude 
at infinity as a free parameter. Note that for the local- 
ized states the normalization has been chosen to be 1 
rather than the number of atoms N because N has been 
absorbed into the coefficient 77. Unlike in linear quan- 
tum mechanics, the nonlinear term in the NLS makes 
the choice of normalization non-trivial. 



A. Localized states 

Localized states are defined as bound states for which 
the tails approach zero at a large distance from the well. 
The functions sech(x) and csch(x) satisfy this criterion 
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and solve the cases of attractive and repulsive nonlinear- 
ity, respectively, i.e. 77 < and 77 > 0. The fundamen- 
tal difference between these two functions and exp(±x), 
which solves the analogous linear Schrodinger problem, 
is not their functional form. Rather it is an offset b which 
for attractive nonlinearity retracts the decaying tail into 
the well and for repulsive nonlinearity pushes the tail into 
the barrier. That is, for x > I the functional form of the 
wavefunction is sech(a; + b) with b < for the former case 
and csch(x + b) with b > in the latter one. This can be 
understood as a physical effect of the nonlinear term in 
Eq. ([j]). Note in the limit as x — > ±00, sech(a; + b) and 
|csch(a; + b)\ both approach exp[=F(x + b)]. 

1. Attractive nonlinearity 



Symmetric solutions take the form 

fx(x) = Asech[fc(x - b)] , (4a) 
f 2 {x) = A 2 cn{k 2 x I to) , (4b) 
h{x) = Asech[k{x + b)\, (4c) 

and antisymmetric solutions take the form 

fi(x) = — Asech[k(x — 6)] , (5a) 
f 2 (x) = A 2 cn[k 2 x + K(m) \ m] , (5b) 
f 3 (x) = Asech[k(x + b)} , (5c) 

where A, k, k 2 , 6, and m are free parameters. fi(x) and 
fs(x) have been chosen in a manner to preserve symme- 
try under the reflection operation x — » —x. In App. it 



is proven that the analogous solutions with dn in place of 
cn in Eq. ([|b) or Eq. (gb) cannot simultaneously satisfy 
normalization and boundary conditions. The dn elliptic 
function is a solution type dependent upon special sym- 
metries, as has been exposited elsewhere [ p0| . The finite 
square well breaks these symmetries. 

Eqs. (^) and Q are subject to the same solution 
method. Consider first Eqs. (||). Substituting them into 
Eq. (Q) with the potential Eq. (||), the conditions 



A 2 = -2k 2 J j] (6a) 

/j = V a -k 2 (6b) 

A\ = -2mkl/r) (6c) 

fi = -(2m- \)kl (6d) 



are obtained. The derivatives of elliptic functions, used 
to obtain Eqs. (^), may be found in the literature p3[ . 
The boundary condition fi(—l) — fi{— I) is equivalent to 
/3O = for Eqs. (H), and requires 

Asech[k(l + b)} = A 2 cn(k 2 l | m) . (7) 

Continuity of the first derivative requires 

Ak tanh[fc(Z + b)] sech[fc(/ + b)] 

= A 2 k 2 <hi(k 2 l j m) sn(k 2 l \ m) . (8) 



Finally normalization requires 

2A\ / dxcn 2 (k 2 x I m) 
Jo 

/oc 
dxsed^ikix + b)} = 1. (9) 

Using the Jacobian elliptic identity m cn 2 (x \ m) + 
(1 — m) = dn 2 (a; | m) and the integral relation 

J Q K ^ dxdn 2 (x I m) = E(m), together with the hy- 
perbolic relations dy sech 2 (y) = 1 — tanh(y) = 
exp(— y) sech(y), Eq. (^|) can be written in the form 

2A 2 

\E(k 2 l I m) - (1 - m)k 2 l] 

k 2 m 
2A 2 

+ — exp[-fc(Z + b)} sech[fc(/ + 6)] = 1 . (10) 

where E(k 2 l \ m) is standard notation for an incomplete 
elliptic integral [p3[ . 

Equation of Eqs. (^b) and (|6]d) and substitution of 
Eqs. (|a) and (§c) into E q s - &, (§)> and © produce 
a system of four simultaneous equations: 

(1 - 2m)k 2 l 2 + k 2 l 2 = V l 2 , (11a) 

klsech[k(l + b)] = y/m k 2 l cn(k 2 l \ to) , (lib) 

k 2 l 2 tanh[fc(/ + b)] sech[fc(Z + b)] 
= ^/m k\l 2 dn(k 2 l \ to) sn(k 2 l \ to) , (He) 

k 2 l[E(k 2 l I to) - (1 - m)k 2 l] 
+kl exp[-k(l + b)]sech[k(l + b)} = -l<q/A. (lid) 

All equations have been multiplied through by an extra 
factor of /. A notational simplification can be made by 
defining a = kl, (3 = k 2 l, u = k(l + b), and 7 = \fV l. 
Then 

(1 - 2m)0 2 + a 2 = 7 2 , (12a) 

asech(w) = v / ro/3cn(/3 | m) , (12b) 

a 2 tanh(w) sech(w) 
= ^fm,fi 2 dn(/3 | to) sn(/3 | to) , (12c) 

(3 [E{j3 I to) - (1 - m)(3] + a sech(w) 

= -lf)/4:. (12d) 

These four equations are in four unknowns, a, (3, u>, and 
to, with an additional three variables I, rj, and 7 which 
are externally determined experimental parameters. 

This system of equations may be compared to the anal- 
ogous set derived for the linear Schrodinger equation, as 
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may be found in any undergraduate quantum mechanics 
textbook [Q. In the linear case, the energy condition 
Eq. (p"l|a) describes a circle. Here it describes a hyper- 
bola for 1 — 2m < and an ellipse for 1 — 2m > 0. Equa- 
tions (|lT|b) and (|TT|c), derived from boundary conditions 
on f(x), are similar in form to the linear ones, which are 
usually combined by dividing one by the other. Finally 
Eq. (|il|d) has no linear analog, as in the linear case nor- 
malization is trivial, and therefore here forms an extra 
condition. The linear system is reduced to two transcen- 
dental equations and solved by graphical methods. Here 
a similiar approach is used. 

Two unknowns may be eliminated as follows. Equa- 
tions dl2^,) and (|l2|b) are solved for a and lu. These are 
then substituted into Eqs. (|T^c) and (p~2|d) to obtain 

y/m tanh{sech _1 [Xcn((3 | m)j} 

= A dn(/3 | to) sc(/3 | m) , (13a) 

v / ro/3cn(/3 | to) exp(— {sech -1 [A cn(/3 | m)]}) 

+/3 [E{/3 | m) - (1 - m)j3] = , (13b) 

where A = \/ra(3/ \/j 2 — (1 — 2m)f3 2 . The resulting two 
equations in the two unknowns (3 and to may be solved 
implicitly by a combination of graphical and numerical 
techniques. The method is demonstrated in Fig. |l|. The 
approximate location of the solution in parameter space 
is first determined by graphical examination of the inter- 
section of Eqs. ([l3|). Then the well-known secant method 
is used to find (to, /?) to double precision accuracy. The 
largest value of to for which Eqs. ( |l3| ) intersect gives 
the ground state and the excitation levels are thereafter 
monotonically decreasing in to. In Fig. [j], a well of depth 
V = 25 and width 21 = 2 and a nonlinearity of rj = — 8 
were used. The intersection in the lower right corner of 
the plot, at to — 0.9258 and (3 = 2.2709, gives the ground 
state. The other illustrated intersection, at m = 0.1934 
and (3 — 4.0749, gives the first symmetric excited state. 

The antisymmetric case Eqs. (|^) may be treated with 
a similiar solution method. The system of equations 

(1 - 2m)(3 2 + a 2 =j 2 , (14a) 

asech(w) = /?cn[/3 + K(m) | m] , (14b) 

a 2 tanh(w) sech(a-i) 
= f3 2 dn[/3 + K(m) \ to)] sn[/3 + K(m) \ to] , (14c) 

/3 {E[/3 + K(m) | m] - (1 - m)P} 

+a sech(w) = -lr)/A (14d) 

results. The ground state and first three excited states, 
two of which are symmetric and two of which are anti- 
symmetric, are plotted in Fig. [2] for r\ = —50, V = 100, 
and 21 = 2. For the ground state, the Jacobi elliptic 



parameter to = 0.999999999972 and the chemical poten- 
tial is /i = —156.25. For the first three excited states, 
to and n are 0.9857, 0.7694, 0.5318 and -40.20, -19.37, 
and -2.982, respectively. Higher excited states are closer 
to linear Schrodingcr equation, sinusoidal-type solutions, 
while lower states are visibly nonlinear. The ground state 
become more bound, i.e. the eigenvalue becomes more 
negative, as the absolute value of the nonlinearity in- 
creases. Atoms tend to cluster at the center of the well 
because of the attractive nonlinearity, as can be seen in 
Fig. @(a). 

2. Repulsive nonlinearity 



Symmetric solutions take the form 

f^x) = -Ac8ch[k(x -b)], (15a) 

f 2 (x) = A 2 sn[k 2 x + K(m) | to] , (15b) 

f 3 {x) = Acsch[k(x + b)}, (15c) 

and antisymmetric solutions take the form 

fi(x) = Acsch[k{x-b)], (16a) 

f 2 (x) = A 2 sn(k 2 x I to) , (16b) 

fo{x) = Acsch[k{x + b)] , (16c) 



where amplitudes and offsets have again been chosen 
to preserve symmetry under reflection. Substitution of 
these equations into Eq. ([!]) with the potential Eq. (||) 
determines the amplitudes and chemical potential 

A 2 = 2k 2 / V , (17a) 
A\ = 2TOfc 2 2 /r ? , (17b) 
fi = (1 + m)k 2 2 = V - k 2 , (17c) 

and, together with boundary conditions, the simultane- 
ous transcendental equations which determine the solu- 
tions up to the three external experimental parameters 
77, I, and 7 are 

(1 + m)l3 2 + a 2 = 7 2 , (18a) 

a csch(cj) = y/m (3 sn[/3 + K(m) \ m] , (18b) 

a 2 coth(w) csch(w) = ~^/m(3 2 dn[[3 + K{m) \ to] 

x cn[j3 + K(m) | to] , (18c) 

(3 {13- E[f3 + K(m) \ to] + E{m)} 

+a e~ u csch(w) = lr//4 , (18d) 

for the symmetric case. A similiar set is obtained for the 
antisymmetric case. 

The ground state and first three excited states are plot- 
ted in Fig. for 77 = 50 and the same well dimensions 
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as were used in Fig. [| V — 100 and 21 = 2. In panels 
(a)-(d) the elliptic parameter m and chemical potential /i 
are 0.9967, 0.8808, 0.6588, 0.4696 and 29.47, 38.23, 49.40, 
63.45, respectively. As the nonlinearity is increased the 
ground state chemical potential increases until it reaches 
the height of the well for rj ~ 210. For 77 > 210 there are 
no normalizable localized bound states in the potential 
well. The solution then becomes partially localized, as 
described in the next section. 



B. Partially localized states 

Partially localized states are defined as states for which 
the tails approach a constant non-zero value at a large 
distance from the well but the density |/(a;)| 2 is still 
higher in the region of the well. They have no analog in 
the linear Schrodinger equation. As stated in Sec. |, the 
analog of scattering states, which continue to oscillate as 
x — > ±00, are not considered. Rather, solutions which 
follow the adiabatic evolution of the localized, bound 
states described above, i.e. those which continue to have 
the same number of nodes as rj is increased, are presented. 

The ground state of the NLS with V(x) = and 77 > 
is a constant, f(x) = ^/Jl. Partially localized states 
can be understood as a deformation of this ground state: 
the wavefunction spills out over the top of the well and 
approaches a constant non-zero value as x — ► ±00. If 
the dispersion term in the NLS is neglected, in the so- 
called Thomas-Fermi limit (?]] , the ground state solution 
is f{x) = yju — V(x) inside the well. If /i < V a the 
wavefunction is zero outside the well; if fi > V then it is 
a non-zero constant. The ground state analytic solution 
presented here is, in some regimes, similiar to that of the 
Thomas- Fermi limit. But it gives the correct, continuous 
connection between the wavefunction inside and outside 
the well and, unlike the Thomas-Fermi limit, describes 
the excited states as well. 

Symmetric solutions take the form 

= -Acoth[k(x-b)}, (19a) 
h{x) = A 2 sn[fc 2 £ + K{m) \ m] , (19b) 
/ 3 (ic) = A<xtih[k{x + b)], (19c) 

and antisymmetric solutions take the form 

h(x) = Acoth[k{x - b)\ (20a) 
f 2 (x) = A 2 sn(k 2 x I to) (20b) 
f 3 (x) = Acoth[k{x + b)} . (20c) 

The amplitudes and chemical potential are 

A 2 = 2k 2 /rj, (21a) 
A\ = 2mk\/r\, (21b) 
/i = (l + m)fc| = V Q + 2k 2 , (21c) 



while the set of simultaneous transcendental equations 
are 

(1 +m)/3 2 - 2a 2 = 7 2 , (22a) 

a coth(aj) = \fm (3 sn[/3 + K(m) | to] , (22b) 

a 2 csch 2 (w) = —y/mf3 2 dn[/3 + K(m) | m] 

x cn[/3 + K(m) \ m] , (22c) 

for the symmetric case. A similiar set is obtained for the 
antisymmetric case. For partially localized states there 
is no normalization condition. Instead A is a free pa- 
rameter. The ground state and three excited states are 
shown in Fig. for rj = 250, the same well dimensions 
as in Figs, g and 0, and k = 1. The parameter m in 
panels (a)-(d) is 0.9999983, 0.9948, 0.9272 and 0.7466, 
respectively. Here k is a free parameter that determines 
the chemical potential and the amplitude A. A is also 
determined by 77, as is apparent in Eq. (|21|e 



C. Transition from localized to partially localized 

states 

As mentioned in Sec. ffl, localized and partially localized 
states are differentiated by the chemical potential. For 
V > 0, [i < V results in the former solution type and 
fi > V in the latter one. It is therefore natural to verify 
that they converge in the limit in which /i — ► V . It is 
once more emphasized that the transition is an adiabatic, 
continuous one which fixes the number of nodes. 

Consider first the limit of the partially localized states. 
(i — > V Q + implies that A — > 0, which in turn implies k — > 
0. Recall that a = kl and u> = k(l + b). This implies that 
the arguments of the hyperbolic functions in Eqs. (p2) 
approach zero. In the limit as y — > 0, coth(?/) — ► ljy 
and csch(y) — > 1/y. Therefore acoth(w) — * 1/(1 + b), 
a 2 csch 2 (w) — > [1/(1 + b)] 2 , and to zeroth order Eqs. ( jgg ) 
become 



(1 + m)l3 2 = 7 2 , (23a) 

l/(l + b)=y/mPsa.[P + K(m)\m], (23b) 

[1/(1 + b)] 2 = -^E/3 2 dn[/3 + K(m) \ m] 

x cn[/3 + K(m) \ m] . (23c) 

Cons ider ne xt the limit of the localized states presented 

v- 



in Sec. II A 2, /i —>■ V~. In this case Eq. (|17|c) requires 
k — > 0. Then by the same reasoning as was used for 
the partially localized state limit above, to zeroth order 



Eqs. become 



(1 + to)/3 2 = 7 2 



(24a) 
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1/(1 + b) = V^/?sn[/3 + K(m) \ m] , 



(24b) 



III. LIMITING CASES 



[1/(1 + b)} 2 = -Vm/3 2 dn[/? + K(m) \ to] 

x cn[/3 + K(rn) \ m] , (24c) 



(3{p- E[(3 + K(m) | to] + E(m)} 
+l/(l + b) = lrj/4. 



(24d) 



The sets of equations (23) and ( p4| ) are identical up to 
the normalization Eq. (24 4), and, as A — > for the par- 
tially localized states in this limit, normalization may be 
imposed. 

The two solution types are also distinct in their non- 
linearity rj, i.e. there is an r]o for which a transition oc- 
curs for a particular set of well parameters and excitation 
level. This provides a simple physical interpretation, as 
is discussed in Sec. 



IV 



Equations ( |24| ) can be solved im- 
plicitly for to and r\ when reduced to two equations of 
the form 

to = — ds[A + K(m) | to] cs[A + K(m) | to] , (25a) 



r?o = f {A - E[X + K(m) | to] + E(m) 



to sn[A + K(m) | to]} , 



(25b) 



where A = j/\/l + m. In practice, to may be as singu- 
larly close to 1 as to = 1 — 10~ 60 , especially for deep wells. 
Equations ( p5|) may be solved for any excitation level by 
first finding the roots of Eq. fl25|a) in descending order 
and then using the resulting values of to to determine 770 
from Eq. (p5|b), with care given to excluding extraneous 
roots. 

To find the transition for the ground state the appro- 
priate simplification is to expand Eqs. (25 1) and (f25|b) 
to first order in toi = 1 — to, solve for toi in the first 
equation, and substitute the result into the second. By 
this method the ground state transition occurs for 
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m = t 



V2 



-V27 



1 + 2e-^T - e- 2 ^ 



2 7 



(26) 



In the limit for which 7 » 1, 
further reduces to 



i.e. for a deep well, rjo 



Vn 



2V2- 



7 — 7e 



(27) 



which to lowest order is simply 770 = 2Vq- For typical ex- 
perimental parameters, as discussed in Sec. the error 
in this approximation is negligible. 

Because determining rjo requires numerical methods 
and/or approximations, /i is a better analytical quantity 
to characterize the dividing line between the localized 
and partially localized states. Then each of the three 
solution types exists in a separate regime: the localized 
states of Sec. 



II A 1 



for ry < 0; the localized states of 
Sec. II A 2 for 77 > and fj, < V Q ; and the partially local- 
ized states for 77 > and fi > V - Note that the case of 
?/ = is that of the linear Schrodinger equation. 



Two limiting cases are presented. In Sec. Ill A it is 
demonstrated that the finite square well potential gives 
the same localized and partially localized solutions as the 
delta function potential 



V(x) = (5(x) 



(28) 



in the limit I — ► and V — > 00 with the area of the 
square well, 2V I, held constant. Eq. ( pS| ) models a fixed 
impurity or scatterer of a size m uch s maller than the 
condensate healing length. In Sec. Ill B it is shown that 
the limit V — > 00 reproduces the full set of statio nary 
states of the NLS under box boundary conditions f^^of . 
Other limits, such as 77 — * or I — > 00, are not worked 
out explicitly but reproduce the correct solutions to those 
cases. 



A. Delta function 

Consider the delta function limit of the localized, 

In order for the 



bound solutions presented in Sec. II A 1 



limit to be tractable, the square well potential of Eq. (|2j) 
must be shifted by V , so that the potential is every- 
where zero except between I and —I, where it is — V , and 
fi = — k 2 outside the well. Let I — ► and V — > 00 with 
2V l = k held constant in Eqs. (|l~T|) . Equations (|ll"|b) 



to ki 



(114) then reduce to 

k sech(fc b) 

k 2 tanh(fc b) sech(fc b) = y/m k 2 (fc 2 • 

k exp(— k b) sech(fc b) = —77/4 , 



(29a) 



(29b) 



(29c) 



where sn(fc2 1 \ to) ~ I has been used. Dividing 
Eq. (f29|b) by Eq. (|29|a), solving for k 2 l, and substitut- 
ing it into Eq. (pTk) with to — > 1 results in 



(30) 



2fctanh(fc6) 



If one makes the identification k = ( then Eq. ( J3C 
together with Eq. (|29|c), correctly describes the localized 
solution to Eq. (Q) with the potential Eq. ( p8[ ) of the form 



Asech[k(x - b)] x < 
Asech[k(x + b)} x>0 



(31) 



for 77 < 0. Equation (pQ) is the delta function boundary 
condition while Eq. (|29|c) is the normalization, as the 
reader may verify, k and b can be solved for explicitly to 
yield 



-21og e (l+4C/?j) 
')+2C 



h — —i — 2 

K - 2 4 



(32a) 



(32b) 



G 



The amplitude and chemical potential are the same, 
namely A 2 = —2k 2 /ri and /i = — k 2 . Note that in the 
limit in which r; — > the linear Schrodinger equation is 
recovered, k — — C/2, and /i = C 2 /4, which reproduces 
the eigenvalue of the well-known linear solution. 

T he delt a function lim it of the solutions presented in 
Sec. [I A 2| and Sec. [II B| is taken in a similiar manner. 
These are similiar to Eq. (|3l|) but have the functional 
form f(x) — ±A csch[fc(x±&)] and f(x) = ±Atanh[fc(a;± 
b)} for localized and partially localized solutions, respec- 
tively. Non-trivial phase solutions, which are essentially 
captured solitons, have been treated by Hakim |27| . Note 
also that Atax\h(kx) is a solution since it has a node at 
the delta function, as well as various Jacobian elliptic 
functions which also have a node at x = 0. These cases 
are not treated here. 



B. Infinite square well 

It sha ll be demonstrated that the localized solutions 
of Sec. II A reproduce the infinite square well solu- 
tions |l2,|20[ in the limit in which V — ► oo. Partially 



localized states have no such limit, as they cannot satisfy 
the boundary conditions. 

Consider first the localized solutions of Sec. II A: 
V — > oo implies that k — > oo. Then Eqs. and 
reduce to 



-2mk 2 /rj . 



-(2m 



= \fm k 2 l cn(k 2 l | rn) , 



(33a) 



(33b) 



(33c) 



k 2 I [E(k 2 l | m) - (1 - m)k 2 l] = -lrj/4 . (33d) 

The same must hold true for I — > — I. Equation (|33|c) 
then requires that the argument of the elliptic function 
cn be an integer number of half periods, i.e. k 2 l = (2p + 
l)K(m), where p is an integer. E[(2p + l)K(m) \ m] — 
(2p + \)E{m) 1 from which it follows that the amplitude, 
chemical potential, and normalization condition are 



A\ = -2m[(2p + l)K(m)) 2 /r] , 



fi = -(2m- l)[(2p + l)K(m)f 



(34a) 



(34b) 



(2p + l) 2 K(m) [E(m) - (1 - m)K(m)] = -\ .(34c) 

In the notation of reference [20], —r\ = 2\~ 2 and jl — 
—/J,/r]. 2j = 2p+ 1, where j is an integer, since the edges 
of the well shift from — I, I to 0,2/. With these identifi- 
cations Eqs. (|34|) correctly reproduce the infinite square 
well solutions for attractive nonlinearity |2^]. For repul- 
sive nonlinearity_[f2l , the limit of the localized solutions 



IV. APPLICATION: TUNABLE TUNNELING 

The application of this work to tunable tunneling fol- 
lows directly from the fact that r\ cx aN. In recent ex- 
perimental developments it has been shown that the scat- 
tering length of certain Bose-condensed atoms, especially 
85 Rb, can be efficiently tuned over several orders of mag- 
nitude, including the sign [pl,p5|]. In Fig. || the effect of 
this tuning is illustrated. For r\ < 0, shown in panel (a), 
the tunneling can be made so small that the wavefunc- 
tion is essentially unaffected by the potential. In panel 
(b) the linear Schrodinger solution is shown, i.e. rj = 0. 
Panel (c) shows what happens when r\ > and \x < V . 
The wavefunction is pushed far out into the barrier, as 
compared to the linear case in panel (b), but is still lo- 
calized. Finally, in panel (d) a partially localized state 
is shown. r\ is increased so much that fj, > V a , and the 
wavefunction spills out over the top of the well and ap- 
proaches a non-zero constant at infinity. A well depth 
of V — 100 has been used for the figure, but the gen- 
eral form of the transition is the same for any well. Note 
that the tuning between panels (b) through (d) do not 
necessarily require a Feshbach resonance, as the density 
can be increased by evaporative cooling, which controls 
the temperature and therefore the number of atoms N 
in the condensate, or by dynamically changing the trap 
potential. 

In Fig. H the localized states shown in panels (a) 
through (c) have a finite number of particles. But as 
the partially localized state shown in panel (d) has a 
constant density as x — > ±oo it requires an infinite num- 
ber of particles. This latter solution type, which has no 
analog in linear quantum mechanics, may be treated ex- 
perimentally in two ways: either the entire system can 
be enclosed by a trap of infinite depth, as for example a 
shallow harmonic potential or an infinite square well po- 
tential with walls far from the finite square well potential, 
in which case it can be tuned by the scattering length or 
the density in the same manner as the localized solutions; 
or it can be controlled by the chemical potential /i > Vo- 
In present experiments the trap depth is effectively infi- 
nite and is usually idealized by a harmonic potential . 
Thus partially localized states have not been observed. 
However, with the advent of BECs trapped in hollow 
blue-detuned laser beams capped on cither end by laser 
light sheets of adjustable intensity it should now be 
possible to observe both localized and partially localized 
states. Although 87 Rb has been the atomic species used 
in hollow laser trap experiments, 85 Rb could equally well 
be used [^6| . Thus tunable tunneling can be studied with 
present BEC technology. 

To explicate the usefulness of this application the phys- 
ical scalings of Eqs. ([!]) and (||) are here worked out 
explicitly. The quasi-onc-dimcnsional Gross-Pitaevskii 
equation with tight transverse confinement and with all 
physical units is 



treated in Sec. II A 2 is found in a similar manner. 
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2m mLyL z 



fo(Xo) 

Hofo(x )- (35) 



This equation is obtained by projection of the three- 
dimensional wavefunction onto the ground state in the 
transverse dimensions, which assumes an approximate 
separation of variables. Details of this a ppr oximation 
are treated elsewhere p^Jl3] , |37|| . In Eq. (|3|), box- like 
transverse confinement at ±L y /2 and ±L z /2 has been 
assumed. Note that the experiment of Ref. |24|] is capa- 
ble of producing these conditions |58]j|] . a is the s-wave 
interaction length which characterizes binary atomic in- 
teractions and N is the total number of atoms. The 
proper scalings between Eq. ( |35| ) and Eqs. (|l]) and (||) 



correlations necessary to describe the decoupling and re- 
coupling of localized condensates. 
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APPENDIX A: PROOF THAT THERE ARE NO 
NODELESS TRIVIAL PHASE SOLUTIONS 

Consider localized, nodeless solutions for attractive 



are then ±i = [h /(2ml 2 )} n, V 
■q = 8TraNl/(L y L z ), x a = lx, and f Q = f/y/l- 

For the case of 10 4 atoms and a trap of dimensions 
21 = lOO/xm and L y , L z — 10/im, a trap of V ~ 10/iK 
depth would correspond to V a ~ 10 8 . The transition be- 
tween localized and partially localized states would oc- 
cur, in this trap, at 770 = 2.0 x 10 8 . Assuming the use 
of a Feshbach resonance rather than a change in den- 
sity to achieve the transition, the scattering length for 
85 Rb, normally ~ —400 a j2j| , need only be tuned to 
a ~ 14 a Q , well within the range of present experimental 
capabilities. 



V. SUMMARY 

It has been demonstrated that the tunneling and lo- 
calization of the macroscopic wavefunction which de- 
scribes the mean field of Bose-Einstein condensates in 
a trap of finite depth is an external experimentally con- 
trollable parameter. Specifically, the full set of localized 
and partially localized stationary states of the nonlinear 
Schrodinger equation with a finite square well potential 
have been presented. Localized states are bound, and 
their tails approach zero at a large distance from the 
well. Partially localized states are defined as states for 
which the tails approach a constant non-zero value at a 
large distance from the well but the density is still higher 
in the region of the well. The latter have no analog in the 
linear Schrodinger equation. Tuning between these dif- 
ferent kinds of states could be studied with present BEC 
technology. 

An obvious extension of these tunable tunneling results 
would be to consider the problem of a condensate in two 
adjacent finite square wells. Analytic solutions of the 
NLS will also be available in this case |3!|. Tuning the 
nonlinearity rj can then increase or decrease the coupling 
between the spatially separated condensates. The NLS 
will not provide a complete description as this coupling 
, but the availability of analytic 
1 provide the zeroth order nonlin- 
ear wave functions needed in models which include the 



[h I (2m I 2 )] V , nonlinearity, i.e. r\ < 0, of the form 



approaches zero J3l| , p2 
mean field solutions wi 



fi(x) — Asech[k(x — b)] , 
f 2 {x) = A 2 dn(k 2 x I to) , 
f 3 (x) = Asech[k(x + b)} , 



(Ala) 
(Alb) 
(Ale) 



where A, k, k 2 , b, and m are free parameters. fi(x) and 
fs(x) have been chosen in a manner to preserve symmetry 
under the reflecti on ope ration x — ► —x. By the methods 
explicated in Sec. II A 1 the two systems of simultaneous 
equations which constrain the parameters are 



A 2 = -2k 2 jx), 



/' 

/' 



V 



-2mfcf/?7, 
-(2 - m)k\ 



and 



-(2-m)f3 2 + a 2 = 7 2 , 

a sech(u>) = j3 dn(/3 | to) , 

a 2 tanh(w) sech(w) 
= to [3 2 cn(/3 I to) sn(/3 | m) . 



(A2a) 
(A2b) 
(A2c) 
(A2d) 



(A3a) 



(A3b) 



(A3c) 



/3E(P I to) 



a e 



'sech(w) = -Z77/4, (A3d) 



where a = kl, (3 = k 2 l, ui = k(l + 6), and 7 = \fVJ,. 

Equations ( A3 ) can be reduced to a set of two simul- 
taneous equations i n the parameters (3 and to, again by 
the methods of Sec. II A 1 . One of these equations, which 
must be satisfied, is 



yj\ 2 -/3 2 dn 2 (/3 I m)dn(/3 | to) 
—to (3 cn(/3 I to) sn(/3 | to) = , 



(A4) 



where A = (2 — to)/3 2 + j 2 . It shall be proven that this 
equation is false. 
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First note that 8 an d 7 are positive definite and < 
to < 1. For Eq. (A4) to be true, cn(/3 | to) sn(/3 | to) 
must be positive, so 



< 8 < K(m) or 2K(m) < 8 < 3K(m). 
Then an absolute value sign may be imposed, 



v/A 2 - B 2 dn 2 (/3 J to) dn(fi | to) 
-to 8 \cia(8 I to)| \sn(8 \ m)\ = . 



(A5) 



(A6) 



Equation (A5) and the fact that < to < 1 imply that 
dn(8 I to) >\cn(8 \m)\. Then 



\/\ 2 - 8 2 dn 2 (8 I to) < to/3 |sn(/3 | m)\ 



(A7) 



Square both sides and apply the elliptic identity to sn 2 (8 \ 
to) + dn 2 (/3 I to) = 1. This implies that A 2 < B 2 . Sub- 
stitute back in A = (2 — to)/? 2 + 7 2 and the condition 
becomes 



(2 - to)/? 2 + 7 2 < B 2 . 



(A8) 



This implies that (2 - to) < 1. But < to < 1. There- 
fore there are no nodeless trivial phase solutions in the 
finite square well for any set of well parameters. Con- 
versely, note that nodeless, dn-type solutions do exist in 
the case of a constant potential with periodic boundary 
conditions |2C[] . 
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FIG. 1. Graphical solutions method: shown are the first 
two intersections of Eqs. (p^|a) (solid line) and (|l3|b) (dashed 
line.) These give the ground state and first symmetric excited 
state, at the lower right and upper left of the plot, respec- 
tively. 




FIG. 2. Shown are the ground state and first three excited 
states for localized solutions with attractive nonlinearity r\. 
The box walls are at x = ±1. 




FIG. 3. Shown are the ground state and first three excited 
states for localized solutions with positive nonlinearity r\. The 
box walls are at x — ±1. 




FIG. 4. Shown are the ground state and first three excited 
states for partially localized solutions, for which the tails of 
the wavefunction approach a non-zero value asu ±oo. The 
box walls are at x = ±1. Here the tails of the wavefunction 
have been chosen to approach 




FIG. 5. Tunable tunneling: shown are solutions for (a) 
negative scattering length, (b) no interactions between atoms, 
(c) positive scattering length and chemical potential less than 
the well depth, and (d) positive scattering length and chemi- 
cal potential greater than the well depth. Thus the tunneling 
and localization of the wavefunction is controlled via external 
experimental parameters. A transition between the localized 
states shown in (a) - (c) and the non-normalizable, partially 
localized state shown in (d) can be characterized by the scat- 
tering length a. For 10 4 atoms of 85 Rb in a trap of dimensions 
100x10x10 fim 3 and depth 10 /iK the transition occurs for 
a ~ 14a . 



TABLE I. Limits of Jacobian elliptic functions and inte- 
grals. The first three, i.e. sn, cn, and dn, are periodic solu- 
tions in the well while the m — 1 limits of dn, ds, and ns, 
i.e. sech, csch, and coth, constitute the decaying tails of the 
wavefunction outside the well. 4K(m) is the periodicity and 
the elliptic integrals K(m) and E(m) both play a role in the 
system of equations which describe the solutions. 

m — m = 1 



sn(u | m) 
cn(u | m) 



sin(u) 
cos(u) 



tanh(ii) 
sech (it) 
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dn(ti | m) 1 sech(w) 

ds(u | m) csc(w) csch(u) 

ns(w | m) csc(u) coth(w) 

K(m) tt/2 oo 

E{m) tt/2 1 
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